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^> Abstract 

(f) , A coarsened model for a binary system with limited miscibility of 

components is proposed; the system is described in terms of structural 
C^D ' states in small parts of the material. The material is assumed to have two 

CO alternative types of crystalline local arrangements associated with two 

components of the alloy. Fluctuating characteristics of a cluster are the 
type and the space orientation of its crystalline arrangement. There are 
two different phase transitions in the model system, an orientation order- 
disorder transition representing melting, and a phase transition between 
phases differing in concentration of components. Depending on the pa- 
rameters characterizing the interaction in the system, this last transition 
q may take place both in the crystalline and in the amorphous (molten) 

, phase. A special approximation is used to study the thermodynamics of 

the system. The calculated phase diagram describes, at least qualitatively, 
the most important features of a binary system. 
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1 Introduction. 



In binary systems, composition-temperature phase diagrams show the regions 
of phase stability, usually under constant pressure conditions pQ. Phase tran- 
sitions in those systems include melting , which is a universal phenomenon in 
all crystalline materials, eventual phase transitions due to non-miscibility of 
components and, polymorphous phase transitions that may appear if there is 
a competition between different crystalline phases of the same material. The 
appearance of these latter transitions depends on the atomic scale interactions 
in the material. In isomorphous alloy systems (e.g. copper -nickel alloy), only a 
single type of crystal structure is observed. A more complicated picture appears 
in alloys with low miscibility of components. In those systems, clusters having 
different predominant components may differ in crystalline arrangements. A 
known example of an eutectic system is the copper-silver alloy, while that of a 
monotectic system is given by the copper-lead alloy. . Similarities of melting in 
different materials allow one to assume that the melting behavior depends on 
rather general characteristics (e.g. the particle size, first moments of the interac- 
tion potential) of the interparticle interactions, and that a coarsened statistical 
model operating with only a few energy parameters may reproduce the main 
features of those materials. In the present paper, we propose such a coarsened 
model in terms of structural states of the system small parts. 

In liquids and solids, the interaction of constituent particles is strong and 
substantially limits the relative positions of particles. The resulting strong cor- 
relation in small volumes of the material is referred to as the local order in 
condensed systems. The term "local" refers to a small part of the material, 
whose size may be estimated from the correlation functions as including one 
or more coordination shells. We refer to this small part as a cluster. To ac- 
count for the limitations imposed on relative positions of cluster particles by 
the local order, one describes the particle arrangements in terms of local order 
parameters. The concept of local order assumes a coarsening of description to 
a space scale larger than the interatomic distance. The resulting model is then 
a simplification that allows one to understand the global ordering in terms of 
the local order. 

In simple liquids (e.g. noble gas liquids), radial correlation functions El 
show some order in the first coordination shell, and a weaker correlation of 
particles positions in the second and higher correlation shells. The possibility 
to describe these liquids in terms of a definite local order is an open question. 
In liquids with stronger interaction of particles, and especially with covalent 
or/and hydrogen bonding, one expects the local order to be stronger, and closer 
to the local arrangement in a crystal. This will be referred to as a crystalline 
local arrangement. 

In the current article, we consider a material in which the local order is 
maintained by a strong short range interaction, so that this local order does not 
significantly change even when the global order is lost by melting. In such a 
material, the loss of the global order during melting results from proliferation 
of topological defects, with the density of these defects small enough to allow 



2 



recognition of the crystalline local structure. In the crystalline state, orienta- 
tions of all clusters almost coincide, while in the liquid the orientations become 
uncorrelated above some finite distance. For a single component material or an 
isomorphous system with only one type of local order, a model of melting based 
on the above physical picture was studied earlier 0]. A more general theory of 
melting that accounts for amorphous inherent structures is studied in a recent 
paper 

In a crystalline two-component alloy with limited miscibility of components, 
one deals at least with two types of local order corresponding to phases with 
the predominance of one component. In the materials considered, the short- 
range interaction energetically favors, for each sort of atoms, the same sort of 
particles in the nearest neighborhood. The physical picture of the alloy is then 
a cluster-scale solution. Close to the phase separation caused by non-miscibility, 
the material is a solution of one type of clusters in the matrix of the other type 
of local structure . This suggestion determines the choice of cluster states in 
our model. 

In what follows, we propose a model for an alloy based on the above physical 
picture. The meta-Hamiltonian of the model is written in terms of a local order 
parameter describing two possible predominant configurations in small elements 
of the material. The particle arrangement in such a small clement is assumed to 
be of crystalline type, with orientations of these crystalline arrangements being 
the fluctuating characteristics. We study the phase transitions and the phase 
diagram in this model. 

2 The local state model of an alloy. 

We divide the material in equal small parts (clusters) of cubic shape. The cen- 
ters of these clusters form a simple cubic lattice with sites having coordinates r. 
In accordance with the discussion in the previous section, we assume that the 
relative positions of cluster particles are, within the accuracy of small thermal 
displacements and possible packing defects, of a crystalline type. The crystalline 
local arrangement is anisotropic. One can characterize the space orientation of 
this anisotropy by introducing local crystalline axes. In the crystal state, the 
orientations of these axes are close to the corresponding orientations of some 
global axes characterizing the crystal on a macroscopic length scale. In a non- 
crystalline state, the orientation of local crystalline axes become, together with 
the type of the local crystalline anisotropy and the local composition, a fluctu- 
ating degree of freedom. In our model, we consider only this degrees of freedom. 
The universal contribution to thermodynamic properties coming from the small 
oscillations of atoms in the vicinity of an energy minimum configuration may 
be found, for a classical system, from a rather general consideration :6:. 

To simplify the modeling, we coarsen the orientation description in a way 
proposed in 0] [5]. Namely, we divide the orientation space of a cluster in 
orientation cells comprising only slightly differing orientations, and characterize 
the orientation state of a cluster by the number i, i = 1, N of the orientation 
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cell. The size of an orientation cell, and thus the number N of orientation 
states, will be discussed later. We assign to a cluster situated at the point r 
|7] an TV-component vector <Jj{r) by the condition that <Jj{r) = 1 when the 
orientation state of the cluster is j, and <Jj{r) = otherwise. The global state 
of the material is characterized by the configuration {<Jj{r)} and described by 
the field (Tj(r) at all points r. Then, a statistical theory of the material may be 
formulated in terms of the field aj (r) of the local order parameter. 

One writes the probability w{<Ji(r)} of a configuration {<7j(r)} as w ~ 
cxp[— H{ai(r)} /ksT] where ks is the Boltzmann constant, and H{<7i(r)} is the 
meta-Hamiltonian of the system; below, we choose units in which fcs=L The 
technique of meta-Hamiltonians (effective Hamiltonian, nonequilibrium ther- 
modynamic potential) in the coarsened statistical description of condensed sys- 
tems was developed in the theory of critical phenomena (see, e.g. HOE] ). To 
obtain H{ai(r)}, one has to take averages, with Gibbs measure, over all de- 
grees of freedom except for those described by <Ji(r). The meta-Hamiltonian 
H{<Ji(r)} thus gives the free energy of the system in a global state characterized 
by given local states <7i(r) of each cluster. The the meta-Hamiltonian H{<Ti(r)} 
depends, through coefficients defining the mathematical form of H{ai(r)}, on 
temperature T, pressure P, and other thermodynamic coordinates of the sys- 
tem. One may use these coefficients as new thermodynamic coordinates instead 
of P and/or T. To be specific, we consider a system at a given pressure P; 
the thermodynamic coordinates of the state are the temperature T and com- 
ponent concentrations. A change of the pressure will change, in general, the 
meta-Hamiltonian. 

The form of the meta-Hamiltonian follows from rather general physical ar- 
guments jHllHl- One expects H to include independent contributions from single 
clusters (one-cluster part), each such contribution on depending only on the 
state i of the corresponding cluster . The interaction of neighboring clusters 
is described by terms depending on states of more than one cluster. Here, we 
only consider two-cluster interactions. The reason is that the interacting poten- 
tials acting between particles have a range shorter then the cluster size, so the 
interaction is strong only for clusters sharing a border. Then, for a system in 
homogeneous external conditions, the general form of the meta-Hamiltonian is 



Here and below, the Einstein's summation rule is supposed for repeating indices. 
Let us define the field uii(r) as 



the partition function Z(ai(r),T) is a sum over all local states in all clusters. 
The quantity Wj(r) defined by (J2J) is the probability that the cluster centered at 



#M = -I>(0AMr-r>i(O -I>(r)«i(r). 



(1) 




<cn(r) >=w 1 (r) = ^£{ CT } Oi(r)exp(— 
^ = E W exp(-^); 



(2) 
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point r is in the state i , it yields the normalization condition 

N 
i=l 

The partition function Z is the generating functional for the field u>i(r): 

<*(0=Ti^. (4) 
oai(r) 

To analyze the thermodynamic behavior of the system, one can make a 
change of variables in the functional integral defining the partition function Z. 
The original meta-Hamiltonian H{ai(r)} depends on a discrete variable <7j(r). 
We use the Hubbard - Stratonovich transformation to re-express the partition 
function in terms of a continuous lattice variable tp. Let us suppose that the 
matrix My-(r) has the form: 

M ij (r)=E ij J(r). (5) 

where Eij is a matrix in the local state space, and J(r) describes the dependence 
of interaction on distance between clusters. Introducing a real, iV-component 
"conjugated" field tpi(r), one rewrites the expression for Z as 

Z = J IU DMr) ex P {-^ E nr < J~ l {r ~ r')E^ t {r)^{r')}x 

Here, the functional integration Dtp goes over all components tpi of the field ipi(r) 
. The integration path in the complex ■0-space may, for a non- positive definite 
matrix Mij(r), deviate from the real axes to make the integral convergent. The 
advantage of the form © is that now the summation over local states can be 
easily taken: 

F=^ rtr ,J- 1 (r-r')E ij tp i (r)i; j (r')- (7) 
-TEMEieMTi E Mr)+ ai (r))}}. 

The new variable ipi(r) is a continuous one, so one can study the above expres- 
sion in the saddle-point approximation, in which the only allowed configuration 
is the saddle point one. In the next approximation one studies then the fluctua- 
tions in the vicinity of the saddle point configuration, to find the first fluctuation 
corrections to the saddle-point approximation. The saddle-point approximation 
is known to fail at critical points; for phase transitions of the first Ehrenfest 
order, the saddle-point approximation may be expected to adequately describe 
the thermodynamics of the system and phase transitions. In the present article 
we shall restrict ourselves with the first saddle-point approximation. 

The saddle-point configuration ip(r) minimizes the thermodynamic potential 
F in 0, it obeys the equation: 

£ r , J-Hr - r')i>i(r') = Z^ 1 exp[i(^^(r) + a,)], 

Z = E i exp[^(^(r)+a i )]. ( °> 
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To obtain the generating functional Z one substitutes the solution of this equa- 
tion in @. For the probability Wi(r) = —SF/Sai(r) one has 

Wi (r) =< <7i(r) >= Er' •'"V " r') < Mr') > ■ (9) 

In the saddle-point approximation, the quantity F in formula J7J) plays the role 
of the nonequilibrium thermodynamic potential [5| ; at equilibrium, F has a 
minimum. We define LUi(r) for nonequilibrium values of ifi by the formula @. 
Then, the nonequilibrium potential F may be written in terms of u>i(r) 

F = | Er,r' M ij( r ~ r')w i (r)w i (r')- 
-T£ r ln[E, exp(i(E r , Afy(r - r>,(/) + a*))]. UUJ 

The quantity u>i(r) may be interpreted as the probability to find the cluster in 
the i-th state. Because of this probabilistic character of the quantity LOi (r) , we 
will call the saddle point approximation the mean-probability approximation 
(MPA). The MPA approximation is given by the equations 

3^F)=0> (11) 

with an additional condition @. In the explicit form, this equations may be 
written as 

^:(r) = Z- 1 (r)exp(^) 
Ei(r) = Er' M*(r - r')^-(r') + o,(r) (12) 
^)=ELexp(«). 
Here, —Ei(r) may be interpreted as the energy of the cluster when the states 
of surrounding clusters are characterized by the mean probability u>i(r'). If the 
field cti is spatially homogeneous (cti(r) = at), the quantities Ei,u>i are, too, 
spatially homogeneous: 

E t = + on, Mij = E r Af tf (r). (13) 

The normalization condition © holds for solutions of (|12l) automatically. The 
equations (|12f) have, in general, more than one solution; the number of solu- 
tions depends on the values of the parameters My, on, T. As mentioned above, 
the characteristics Mij , on depend on thermodynamic coordinates of the state 
(temperature, pressure etc.), and may be used as new thermodynamic vari- 
ables. Different solutions correspond to different phases of the system. In the 
vicinity of the first order phase transition, one expects at least two minima of 
F, corresponding to the stable and metastable phases, with the lowest minima 
corresponding to the stable phase. By specifying the form of the matrix My, 
one arrives at different models. Some of these models are discussed in following 
sections. 



3 Melting of a single-component material, and 
the Potts model. 

In this section, we discuss a coarsened model for melting in a single-component 
material E] • The discussion serves as an introduction to the modeling in the 
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next section; a more detailed study of one-component melting may be found in 
ref.|S]- Here, the material is supposed to preserve the same crystalline local order 
both in the crystal and in the liquid state. The only degrees of freedom described 
by the model are the local orientations. To give a coarsened description of 
orientation order-disorder, one divides the cluster orientation space in angular 
cells. By definition, all orientations within a given angular cell correspond to the 
same cluster orientation state, labelled by the angular cell number. The number 
of cluster orientation states is then determined by the size of an angular cell. 
With increase of the angular cell size, the description becomes more coarse. We 
choose this size to model the main features of the system in a most schematic 
way. 

In accordance with the experimental data, the ground state ( the state hav- 
ing the lowest energy) is a crystal, with all clusters occupying the same angular 
cell. The globally amorphous (but locally crystalline in materials under con- 
sideration) global state of the material has a higher energy. This implies an 
angular attraction between neighboring clusters: the meta-Hamiltonian is at 
minimum when their orientations coincide and all clusters form an ideal crys- 
tal. One could expect that this angular attraction has an angular size, so that 
the two-cluster contribution to the meta-Hamiltonian increases sharply with the 
disorientation of the clusters, and then substantially flattens when the mutual 
disorientation becomes larger than the angular attraction size. We choose the 
size of the angular cell equal the angular attraction region. Then, a simple 
schematic way to model the energetics of neighboring cluster interactions is to 
assign an angular binding energy J < to two neighboring clusters belonging to 
the same orientation cell (having the same orientation state). The interaction 
energy of two neighboring clusters having different orientation states is chosen 
as the origin of the energy scale E = 0. The resulting model with only one 
angular energy characteristics has the form of an N state Potts model with the 
matrix My having the form 



The assumed equivalence of all orientation states leads to a high symmetry . The 
parameter m is now state-independent and may be put equal to zero without 
loss of generality. In the symmetric high-temperature phase, all N equivalent 
states have equal probabilities 



At lower temperatures, this symmetry is spontaneously broken. At these tem- 
peratures, one of the states, which we label by i = 1, is more probable than the 
others ( uj\ > 1/N), while N-l states i=2,... ,N have equal probabilities 



My = S tj J(r) , 



(14) 



Ui,h = -?f, i = l,-,N. 



(15) 




1 — u>x 



(16) 



N - 1 
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From the minimum condition for the non-equilibrium potential F (see (|1U|) ^ 
one obtains the MPA equation |Q (see ) 

wi = Z- 1 e JvUl ' T , Z = e J ^/ T + (JV — i) e J "(i-« 0/("-i)t (17 ) 

Here, z^J(O) = ^2 r J(r), with z/ being the effective number of cluster neigh- 
bors. Depending on the temperature, the MPA equation has different solu- 
tions describing the high-temperature and the low temperature phases. In the 
high-temperature phase, all states have equal probabilities uji.h — 1/N, so that 
the quantity x = uj\ — jj equals zero. This solution [x — 0) formally exists 
at all temperatures but corresponds to a minimum of F, and thus to a sta- 
ble or metastable state, only at T > T 2 . The low - temperature solution 
x w (AT — 1)/N exists at T < T\, The transition from the low-temperature 
(x 7^ 0) to the high temperature, symmetric phase (x = 0) takes place at 
T = T m , T\ > T m > T2] this transition is of the first Ehrenfest order (discon- 
tinuous). For N 3> 1 the temperatures T±,T2,T m may be found in explicit 
form: 

Tl ~ W T ' 2 ~ Tm ~ 2AlnA ■ (18) 

The model may be applied to chemical compositions of stoichiometric com- 
pounds behaving like pure materials |Hj • The necessary condition of applicabil- 
ity is that the local order should not be destroyed via melting, including the 
chemical order (the material should not dissociate). 

In multicomponent materials there may be a competition of several types of 
local structure. In a binary system A — B , these competing structures may 
correspond to pure components and to all intermediate stoichiometric compo- 
sitions. The existence of competing structures enhances the set of local states. 
In the following section , we generalize the above model to account for two 
competing types of local order. 



4 System with two competing local structures. 

In this section, we propose and discuss an idealized model of a material with 
two possible local structures. To account for the type of local structure as a 
new degree of freedom, we introduce two indices for labelling the local state. 
The vector a is now written as 

a{; j = 1,2. (19) 

The upper index corresponds to the type of local order, while the lower one 
numerates the orientations of a cluster. As in the previous section, we assume 
division of the orientation space into orientation cells, with n orientation states 
for the first and m for the second local structure. The meta-Hamiltonian of the 
model may be written in the form 

n 

ff = - Q EE^M-E^M M ^ r - ; W( ; )- ( 2 °) 

r i—1 r.r' 
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The parameter a describes the internal energy difference between two local 
structures. As in the case of a one-component system, we choose the angular 
cell sizes to comprise the angular attraction regions. The kernel Mj(r) of 
interaction is assumed to differ from zero only for nearest neighbors, and to 
have the form 

M]l (r) = J l (r)5 a ; Mf (r) = J 2 {r)5 a ; M% (r) = Mf (r) = i(r) (21) 

The interaction between clusters of same composition is chosen in the Potts 
form discussed in the previous section. The interaction between clusters having 
different local structures (the term M 12 and M 21 ) is assumed to be orientation- 
independent; the argument here is that if clusters with different structures share 
a boundary, then the energy of the boundary s is mostly due to the difference in 
the neighboring structures and less due to their relative orientations. Parameters 
Ji, J 2 characterize then the orientational interaction of same type clusters, and 
e is an inter-structural "surface" energy. We will denote 

e = je(r)dV; Ji = J Jx{r)dV\ J 2 = J J 2 (r)dV. (22) 

^From symmetry arguments, one has for the mean values of the parameter 
u> (see the previous section): 

Here, x is the concentration of clusters having the first structure. The quanti- 
ties wiand <x>2 in l|23|l describe the orientational ordering; the lower index 1 labels 
the preferred orientation in the orientationally ordered (crystalline) state. Now 
one is able to classify the possible phases in the model: 

1. uji ~ x, x ~ 1 - crystal with the first type of local structure (CI); 

2. lo\ — x/n, x ~ 1 - liquid with the first type of local structure (LI); 

3. 0J2 ~ 1 ™ x, x ~ - crystal with the second type of local structure (C2); 

4. LU2 — (1 — x)/m, x ~ - liquid with second type of local structure (L2) 

The nonequilibrium thermodynamic potential F defined by formula (|10|l may 
be calculated using the meta-Hamiltonian of the model; some lengthy algebra 
then leads to the form 

t - -giyi^i + Ji—^i — 1- ^2^2 + J 2 — ^1 — )+ 

+ex(l -x)- Tln{exp(i Wl + f(l - x) + f )+ 

+(n - 1) exp(4^il + f (1 -x) + f )+ (24) 
+ exp(f uj 2 + fx) + (m - 1) exp( M T \- m x ^ 2) + fx)} 

The minima of this function determine equilibrium phases and the phase dia- 
gram of the model; we study the phase boundaries in the following section. 
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5 Binary system with limited miscibility. 



Consider a binary system in which the two local structures correspond to phases 
with predominance of component A and B respectively. If these phases are sto- 
ichiometric compositions, the model describes the material in some range of 
intermediate concentrations; we assume that no other stoichiometric composi- 
tion can add a third type of local order. We consider systems in which the two 
types of local order differ sharply (have different groups of local symmetry or 
incompatible interatomic distances). Systems with limited miscibility usually 
exhibit a sharp decrease of melting temperature with deviation from the stoi- 
chiometric compound ll(Jj . The experimentally investigated phase equilibria 
and the phase diagrams in these materials may be classified as follows |121 IT5| : 

1. Systems with limited miscibility both in the liquid and solid state (the dia- 
gram with eutectic and monotectic equilibria), for example Pb — Zn, Cu- 
rb, Ga-Pb. 

2. Systems with limited miscibility in the solid state and unlimited miscibil- 
ity in the liquid one (the diagram with single eutectic equilibrium), for 
example Ni — Cr, Cu — Ag, Cu — CuMgi. 

The concentration of the second component determines the volume parts x 
and 1 — x of clusters having the 1 and the 2 structures respectively. Let us study 
the phase diagram of the model in the (x, T) plane (Fig.l). 

The MPA equations may be written as 

wi = f = {1 + (n - 1) exp(Axw 1 i^f )}" 1 
£2 = & = {1 + (m - 1) exp(f (1 - x)!^)}- 1 (25) 
x = i 1 + 37 ex P( zK 1 " x )£>2 - t-^i ~ f (! - 2x ) - 

When x — or x = 1 , these equations coincide with those of the Potts 
model for the melting of pure components. The apparent effect of the second 
component is the lowering of the melting temperature. Besides that, the redis- 
tribution of components between liquid and solid phases [B] results in two-phase 
domains (L2, C2) and (LI, CI) which are restricted by the lines x^{T), x^(T) 
and X5(T),xq(T) respectively (see Fig.l). The boundaries of a two-phase do- 
main may be found by standard procedure from equilibrium free energies of both 
coexisting phases; to find the equilibrium free energy F eq (x,T), one applies the 
method of Lagrange multipliers: 

-F eq (x, T) = J 2 (l - xf{u 2 -|- ^gf) + Jrx 2 ^ - f - |#^)+ 
+sx(l — x) — T{1 — x) ln[aj2(l — x)] — Txln[ajix]. 



The phase equilibrium conditions for, e.g., the domain (L2,C2) are 

p (T rr \ 9F cq C2 (x 3 .T) _ p ( rp \ dF sq} 1,2(3:4, T) . 
feq,C2(l , £3 j X 3 — £ e q,L2\-L , X<l) 9^ X4, 

9F eq ,c2(x 3 ,T) _ dF cqiL2 (x 4 ,T) 



(26) 
(27) 
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Here, the indexes C2 and L2 refer to the phase. For the case to, n 3> 1 , these 
conditions may be obtained in explicit form: 



f (1 " 4) + exl - ±x\ - Tln(l - x 3 ) = 
&(1 - x\) + - ^ - Tln(l - x 4 ) + Tlnm ; 



Tin + ^(1 - Xa) + 2x 4 e - & 



■ fr 4 



At x = or x — 1 these relations give (I18[) for the melting temperature of pure 
material. For a small miscibility [x <C 1) one has approximately 

W =h, —=k Al (29) 

with (to, n» 1): 

2 In to 



The conditions determining x§ (T) and xg (T) may be found in the same way. 

The miscibility gaps both in liquid and solid states and the occurrence of 
two-phase domains (Ll,L2)and (C1,C2) are due to the "interface tension" 
e between two structures. The width of the domain (L1,L2) is limited from 
above by the gap (fig. 2). From the above equations one obtains the relations 
between the parameters of the gap T, x, x* (see fig. 2) as 

F(x) = F(x*), 

F ( x ) = -it* 1 + t x2 + £x2 - Tln (i - x ) ; (31) 

T In £^±f = ( A + i _ 2e)(x* - x). 

(1— x*)x V m n I V / 

Above the domain (L1,L2) in the phase plane the system is homogeneous. 
The line separating this domain from the homogeneous phase is a line of phase 
transitions. The phase transition is of the first Ehrenfest order anywhere except 
the top of the gap with coordinates x c , T c in the x — T plane 

x c = -; T C =A + A_£ (32) 
c 2' c Am An 2' V ' 

where it is of second order. The two-phase domain (CI, C2) is limited by 
lines X\,X2 (see fig. 3). At low temperatures, or for large to, n 3> 1 one may 
approximate the equations as 

-(! - X D - £A + £x i -T\n(l - xx) = 



2 ^ "1J 2n 

-±x% + ^(1 - xl) + ex\ - Tin 



Ml-x l )-^x 1 + 2ex 1 -T\n^^l= ^ 
= &(l-x 2 ) - Jix 2 + 2ex 2 - Tln±^. 

Continuing all lines obtained until their crossing points and using the Gibbs 
phase rule, one gets the full phase diagram of the system. Fig. 4 schematically 
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shows the phase diagram of a system with monotectic and eutectic equilibriums. 
If the crossing point of lines X4, xq lies in the homogeneous phase, the phase 
diagram degenerates into that with single eutectic equilibrium (see Fig. 5). The 
values of x a ,Xb, x*,T* (see fig. 5) are related by 



F eq M^T*)- dF ^- T ^ x a 



- F r (<r* T*\ — 9F ei,h(x* ,T") * _ 
— r eq,L\.x ,±) dx , x — , . 

dF sq , C 2{x a ,T) _ dF^, L {x*,T*) _ dF eq , C i(x b ,T') 
dx a dx* dx}j ' 

where the thermodynamic potentials F eq for corresponding phases are defined 
in (HIJ). 



6 Polymorphous Materials 

The model may be also used, with minor changes, to describe polymorphic tran- 
sitions in a single-component material. In this case two different local structures 
are the competing local arrangements of the same atoms, so the parameter x is 
no longer conserved. The thermodynamic potential (|24[1 and the MPA approxi- 
mation i|25[l give the phase diagram in the (a,T) thermodynamic plane (Fig. 6). 
Starting from one of the phases one changes the parameter a (the difference of 
local energies of the two arrangements playing the role of a thermodynamic co- 
ordinate of the state), and arrives on a polymorphic phase transition C2 — > CI 
at a w J2 — J\ in the solid state. This line may have a continuation in the 
liquid area, where it terminates in the critical point: 

oi c = A -e + T c ln^ -2T C , or . 
T C = I[A + A_2 £ ] 

c 4 L n m 1 

Below T c , the local orders are well distinguished, and phases LI and L2 are 
divided by the line of the first order phase transitions with a critical point at 
the end of the line (see the line AB in Fig. 6). The polymorphic liquid- liquid 
phase transitions were predicted and theoretically studied earlier recently, 
this prediction was confirmed by experiments. |22 ) [T§ | The discovery of liquid- 
liquid phase transitions in single-component materials gives a strong argument 
in favor of a local order in at least some melts. 



7 Discussion of the theory 

The model proposed in this paper is rather schematic, it operates with only a 
few parameters characterizing the material, and does not pretend to give a pre- 
cise quantitative description. At the same time, it gives a qualitatively correct 
description of a surprisingly wide range of known phenomena, and in particular 
predicts the topology of the binary system phase diagram. The proposed model 
may be compared with other known models for the melting transition. 



12 



An early model of melting is the Landau theory with a scalar order parame- 
ter (0, see also in which a cubic term is postulated to result in a first order 
phase transition. In this approach, the nature of the transition and the short- 
range order in the material are not considered. On the opposite extreme, one 
can calculate the thermodynamic potentials for the liquid and crystalline phases 
using "realistic" interatomic interaction potentials , by numerical techniques. 
Our approach is intermediate. The local order in the material has to be found 
from molecular-scale calculations (or from experiment). This local order is sup- 
posed to be a feature both of the crystalline phase and of the melt. The pos- 
sibility of crystalline local order in the melt is not expected to be justified for 
all liquids; the conditions for such a feature is discussed in Both strong 
covalency and local potentials of three (or more) body type (steric interactions, 
hydrogen bonding, ...) favor the local crystalline local order in the melt.. For 
materials with crystalline local order both in the solid and in the molten states, 
the phases may be described in terms of the local structure, using the orientation 
order parameter. This also allows one to address the problem of polymorphic 
phase transitions in the melt (see sec. 4. 2). 

A special situation appears when the competing structures have substan- 
tially different densities. The change a of the internal energy upon changing 
the type of the local structure then includes the large work against the external 
pressure. In such materials, melting may coincide with a polymorphic transition. 
This may explain the "anomalous" melting of covalent crystals with diamond 
structure (C, Si, Ge) and compositions of the type AjuBy (InSb, GaSb) with 
sphalerite structure. The melting of these substances at modest pressures is 
accompanied by a change of the local order, and the coordination number has 
nearly 100% increase 

[El EH . The large density increase causes the semiconduc- 
tor - metal transition during melting 16 . At high pressures, these substances 
transform into metallic crystals with the structure of white Sn |18j . One expects 
that there exists a line of the semiconductor - metal transition which crosses the 
melting line of crystal melting. At this point there is a three - phase equilibrium 
of the covalent crystal, the metallic crystal, and the metallic liquid correspond- 
ing to the point A on fig. 4. In principle, the coexistence line of the crystal and 
liquid having different local structures may terminate (at a lower pressure) at 
the crossing point of the melting line and the line of metallic liquid - covalent 
liquid transition. This line of polymorphic liquid -liquid transition extends then 
into the liquid area and terminates in the liquid -liquid critical point. Recently, 
a liquid- liquid phase transition line of that sort was discovered in selenium |19| . 
A liquid-liquid phase transition with loss of the metallic conductivity was long 
ago predicted by Landau and Zel'dovich|2*U|. 

The "covalent" melting takes place at high temperatures [2JJ- At lower 
temperatures, polymorphic phase transitions in the melt may be expected in 
materials having low-coordination covalent modifications (with the coordina- 
tion number less than 4). Melting of such materials does not require large 
temperatures because the energetics is determined by the non-covalent bonds. 
Then, at high pressures, the covalent melt must transform into a metal one. 
Recently, liquid- liquid transitions were found in liquid Te, Bi and Se |22| . 
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Qualitatively, the phase diagram for Se has the shape shown in Fig. 6. The 
many liquid-liquid phase transitions in Te and Be evidence that there are sev- 
eral competing local structures in the single component liquids. Compositions 
of the type AjyBvi, {Av)2{Bvi)z with covalent melting ^1 (GeTe, SnTe / 
PbTe, PbSe, PbS, Bi^Se^, Bi2Te^, Sb2Te^), may be candidates for liquid-liquid 
metal- covalent transitions. If one supposes that at high pressure any substance 
turns out to be a metal [201 > then this phase transition should take place for the 
elements of V and VI groups, but it is unclear whether it is accompanied by a 
polymorphic transition. The local state representation allows one to model the 
changes of the global order in condensed materials having a "hard" local order. 
This order may have components other than the structure (the local packing 
of atoms, and the orientation of the local structural anisotropy). Similar mod- 
els are being used to describe magnetic [221 IS] an d ferroelectric ordering. 
These general behaviors can, as we have shown, be understood in terms of a 
very simple Potts-type model that includes only a small number of physical 
parameters characterizing the energy difference between local crystalline states, 
and the energy due to intcrfacial interaction between two neighboring clusters 
of same or different types. These parameters are physically meaningful, and 
they seem to be adequate to discuss the phase diagram of binary systems. 



8 Acknowledgments 

We thank Boris I. Shumilo, Antoni C.Mitus and Michael V.Chertkov for many 
discussions of crystalline ordering. This work was supported by NASA (Grant 
NAG3-1932), the Chemistry Division of the ONR, and by the CAMP collabo- 
rative MVRI program of the ONR and by the NSF-DMR through the North- 
western University MRC (Grant DMR 9120521). 



14 



References 

[1] R.E.Reed-Hill,P/i?/szca/ Metallurgy Principles, PWS-KENT Publishing 
Comp., Boston, 1973. 

[2] J.-P.Hanscn, I.R. McDonald. The theory of Simple Liquids, 2-nd edition, 
Academic Press, London 1990.. 

[3] A.R.Ubbelohde. The Molten State of Matter. Wiley, New York, 1978. 

[4] A.C.Mitus, A.Z.Patashinski, Sov.Phys. JETP, 53, 798 (1981); Phys. Lett. 
87A, 179 (1982). 

[5] A.Z.Patashinski, M.A.Ratncr. Inherent amorphous structures and statisti- 
cal mechanics of melting, Submitted to J.Chem. Phys. (May 1996) 

[6] L.D. Landau, E.M.Lifshitz. Statistical Physics, Part 1, Pergamon, Oxford 
1986 

[7] A.Z.Patashinski, M.V.Chcrtkov. Local State Representation in Statistical 
Mechanics of Condensed matter. Preprint INP 91-51, Novosibirsk 1991. 

[8] A.Z.Patashinski, V.L.Pokrovski. Fluctuation theory of phase transitions. 
Pergamon, Oxford, 1979. 

[9] Intermetallic Compounds. Edited by J.H.Westbrook, John Wiley & Sons, 
Inc., New York - London - Sydney, 1975. 

[10] H.J. Goldsmith. Interstitial Alloys. Butter Vorth, London, 1967. 

[11] A.C.Mitus, A.Z.Patashinski, B.I.Shumilo,Phys. Lett. 113A,41 (1985); 
A.C.Mitus, A.Z.Patashinski. Acta Phys. Pol., 38, 10105 (1988). 

[12] R.W.Cahn, P.Hanscn. Physical Metallurgy. North - Holland, Amsterdam - 
Oxford - New York - Tokyo, 1983 

[13] M.Hansen. Constitution of binary alloys. Mc-Graw Hill, New York, 1958. 

[14] T.Kitamura, Physica A, 160, 181 (1989). 

[15] H.O.Carmcsin, Phys. Rev. B, 41, 4349 (1990). 

[16] V.M.Glazov, S.N.Chizhevskaya, N.N.Glagoleva, Zhidkie poluprovodniki. 
Nauka, Moscow 1967 (in Russian); 

[17] A.F.Ioffe, A.R.Regel. Progr. Semicond. 4, 237 (1960). 

[18] G.H. Drickamer. in: Phase stability in metals and alloys, ed. by 
P.S.Rudman, Mc-Graw Hill, New York, 1967 

[19] V.V.Brazhkin, R.N.Voloshin, S.V.Popova, A.G.Umnov, Phys. Lett.A 154, 
413 (1991) 



15 




Figure 1: The thermodynamic (x — T) plane of a binary system. The two - 
phase domains of liquid and crystal coexistence are drawn. 
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Figure 2: The miscibility gap in the liquid state. The position of the gap 
corresponds to eqn. (31). 
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Figure 3: The two-phase domain (Ci and C2) in the thermodynamic plane of 
the system. 
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Figure 6: The phase diagram of a polymorphous system. 
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